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^ : Abstract 

^ ' A theory is presented which quantitatively accounts for the cooperative adsorption of cationic 
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surfactants to anionic polyelectrolytes. For high salt concentration we find that the critical adsorp- 
tion concentration (CAC) is a bilinear function of the polyion monomer and salt concentrations, 
with the coefficients dependent only on the type of surfactant used. The results presented in the 
paper might be useful for designing more efficient gene delivery systems. 
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I. INTRODUCTION 



Polyelectrolyte solutions are ubiquitous. Over the years they have attracted attention of 
Chemists, Physicists, and Biologists. These complex systems have also found many indus- 
trial applications ranging from water treatment to superabsorbants. Besides the charged 
polymers (polyions), a polyelectrolyte solution can contain salt and other molecules which 
may strongly interact with the polyions resulting in a plethora of interesting behavior. The 
complexity of polyelectrolyte solutions and the long range of the Coulomb interaction makes 
study of these systems particularly challenging l^±k,^di^,1^2JJ^dld^dldkd^dldM2i22i2LM., 

In this work we will explore polyelectrolyte solutions containing salt and ionic surfactants. 
These systems have attracted particular attention because the polyelectrolyte-surfactant, or 
more realistically, polyelectrolyte-lipid adsorption can be used to neutralize the DNA charge, 
facilitating the transfection of the DNA across a phospholipid cell membrane. Since the 
polyions which are of particular interest to us, such as the DNA, are quite rigid (persistence 
length 50nm) one can model them as rigid cylinders with the monomeric charge uniformly 
distributed along the main axis. This provides a significant simplification, since for such 
rigid molecules the internal degrees of freedom can be ignored. 

Association between the cationic surfactants and the anionic polyions is driven by the 
electrostatic and the hydrophobic interactions. The hydrophobicity of the surfactant tails 
is responsible for the cooperative nature of the surfactant adsorptionS^i^ii^i. When the con- 
centration of surfactant inside the solution is low, most of the molecules remain free, unas- 
sociated. As the concentration of surfactant is increased, a density is reached at which the 
number of adsorbed surfactants rises sharply. We will denote this as the critical adsorp- 
tion concentration (CAC). The goal of the theory is to predict the value of the CAC as a 
function of the polyelectrolyte and the salt concentrations for different kinds of polyion and 
surfactant molecules. 

All the information about the thermodynamics of the system is contained in its free 
energy. The Helmholtz free energy for the polyelectrolyte-surfactant solution can be con- 
structed based on the fundamental ideas of Debye, Hiickel and Bjerrun>2^i2Si'2Ii2^. The re- 
sulting Debye-Hiickel-Bjerrum theory corrects for the linearization of the Boltzmann factor 
by explicitly introducing into the theory the non-linear configurations in the form of clus- 
ters composed of associated anionic and cationic entities. The concentration of clusters is 
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FIG. 1: The model. 



determined by the law of mass action. It can be shown that in the infinite dilution limit 
this theory yields the exact Manning limiting laws for rigid polyelectrolyte solutions^S. 

To model the hydrophobic interaction of surfactant tails, we will introduce a hydropho- 
bicity parameter x- The basic theory has been presented elsewhere^^i^, and in this paper 
we will review it only briefly, referring the interested reader to the original publications. The 
goal of the present work is to compare the predictions of the theory with the experimental 
data on the association of undecyl-, dodecyl, tridecyl-, and tetradeculpyridinium cations 
with the anionic polyelectrolyte sodium dextran sulfate in an aqueous solutions^. We will 
also derive a simple asymptotic law for the behavior of the CAC at high salt concentrations, 
which we expect to be valid for all rigid polyelectrolyte-surfactant systems. Although dex- 
tran sulfate is not nearly as rigid as the DNA for which the original theory was developed, 
nevertheless, if we restrict ourselves to sufficiently large salt concentration for which the 
Debye length is shorter than the polyelectrolyte persistence length, the polyion flexibility 
can be neglected. We should also note that the present theory is only valid for sufficiently 
dilute polyelectrolyte solutions prior to a phase separation which can occur following addi- 
tion of cationic surfactants. The high density phases will have to be described by a different 
theory^^. 

In the next section we will briefly review the thermodynamic model used to perform the 
calculations. 
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II. THE THEORY 



The polyelectrolyte solution, illustrated in Fig. 1, consists of Np cylindrical anionic 
polyions of length L and diameter Up, monovalent amphiphiles, and Ng "molecules" of 
salt. In an aqueous solution, the anionic polyions are ionized resulting in Z monomers, 
each of charge —q, distributed uniformly with separation b = L/Z along the main axis. 
As a result of ionization, Z monovalent counterions arc liberated into the solution for each 
polyion. We will suppose that salt is a strong electrolyte which is fully dissociated in an 
aqueous environment resulting in an equal number of positively and negatively charged 
monovalent ions. Similarly, we will assume that the head group of each surfactant molecule 
is fully dissociated, producing a negative monovalent colon (charge —q) and a polymeric 
chain with a cationic head group of charge +q. For simplicity, all the counterions and colons 
will be treated as identical, independent of the molecules from which they were derived. The 
ions will be idealized as hard spheres of diameter Uc and charge ±q located at the center. 
The surfactant molecules will be modeled as polymers of Za spherical monomers (each of 
a diameter Uc) with the head monomer carrying a charge +q. The solvent, water, will be 
treated as a continuum of dielectric constant D. The hydrophobic interaction between the 
surfactant tails is short ranged and will be characterized by a parameter x- 

Inside the polyelectrolyte solution, the concentration of polyions is pp = Np/V, the con- 
centration of monovalent salt is Cg = Ng/V, and the concentration of monovalent amphiphile 
(surfactants) is Ca — Na/V. The strong electrostatic interactions between the polyions the 
microions and the cationic surfactants leads to formation of complexes which — in thermo- 
dynamic equilibrium — will be made of one polyion, ric associated monovalent counterions, 
and Ua associated amphiphiles. Particle conservation requires that 

Pc^ (Z - He) Pp + Cs , 

Pa = Ca-UaPp , (l) 

where pc is the concentration of free monovalent counterions and pa is the concentration 
of free amphiphiles. The concentration of colon p_ inside the solution is unaffected by the 
association, 

p_ = + a , (2) 
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The total concentration of free monovalent charges is then 



Pf = Pc + Pa + P- 



(3) 



The goal of the theory is to determine the number of condensed counterions Uc and the 
condensed surfactants n^. In order to achieve this, the Helmholtz free energy of the poly- 
electrolyte solution is constructed and minimized. The details of the calculation can be 
found elsewhere^i^ here we just give the few main steps. 

The largest contributions to the Helmholtz free energy of a dilute polyelectrolyte solution 
are electrostatic and entropic, 

F = F,i + Fent . (4) 

The electrostatic free energy density fei = F^i/V is the result of the interaction between the 
complexes and the free charges inside the solution. It can be obtained using the framework 
of the Debye-Hiickel-Bjerrum (DHBj) theory-!^ and has been calculated in our previous 
work to h^-^: 

el _ PpZUa/L)L, . (kq) 



T*{Kaf 

where (3 = l/ksT, f = F/Np and 

Jo(Ka) 



2 In [KaKi{Ka)] — Xo(Ka) 



(5) 



(6) 



,0 Kf (x) 

Here Kn{x) is the modified Bessel function of order n, n in {naf' = Att p*^/T* is the inverse 
of the Debye screening length, = pfO? is the reduced density, and T* = D ksT a/ 
is the reduced temperature. We have defined a = (oc + ap)/2 as the effective radius and 

the effective charge of the polyion-amphiphile complex. The entropic 
free energy density is given by— 



Pf 
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(7) 



where i G {p, c, a, — } represents the different species and Q is the internal partition of the 
specie i. For structureless particles, the internal partition function is C- = Cc = 1- Foi' 
surfactants, the Flory theory gives Ca+ = ^a- The respective volume fractions (pi are: 
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Here we have introduced the fractions of the polyion monomers associated with the con- 
densed counterions rric = ric/Z and with the condensed surfactants rria = ria/Z. The internal 
partition function of a complex (p can be calculated by modeling the polyion as a one di- 
mensional lattice of Z adsorption sites. The number of associated counterions/surfactants 
at each site can be either zero or one. The problem of calculating the internal partition 
function of the polyion-amphiphile complex then reduces to finding the free energy of a one 
dimensional lattice gas of three different states: empty, associated with a counterion, and 
associated with a surfactant. This free energy can be calculated approximately^ to be 



ln(p [mc,ma\ = 



Z.' 



+ l3x{Z-l) m; 



Z2 

+ Z rric In rric + Z iria In 

+ Z {l-iTLc- ma) In (1 - mc - m^) , (9) 

where ^ = (3 / D a is the Manning parameter, K = Zy^p (Z) — ip (1)] — Z + 1, and ip (n) is 
the digamma function. The first term of Eq. (Q accounts for the electrostatic interaction 
between the polyion and the condensed counterions and surfactants, the second term is due 
to the gain in the hydrophobic free energy when two surfactant tails are in a vicinity of each 
other, the other terms are entropic - resulting from the thermal diffusion of the condensed 
particles along the polyion chain. The characteristic energy of the interaction between the 
condensed surfactants is measured by a phenomenological parameter x- For each type of 
surfactant molecule, the value of x "will be determined by the best fit to the experimental 
data. 

The equilibrium configuration of the system is found by minimizing the total Helmholtz 
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TABLE I: Relation between the fitted x and Za- 
free energy with respect to the number of associated counterions and surfactants, 



Solving this system of equations, the number of condensed counterions and surfactants can 
be determined as a function of the concentration of polyelectrolyte, salt, and surfactant. 



III. RESULTS 



We will define a surfoplex as a complex in which almost all the polyion monomers are 
neutralized by the cationic amphiphiles. Surfoplexes are formed when the concentration of 
the amphiphile inside the solution reaches the CAC, C^. For practical applications in gene 
therapy, the CAC is of particular importance since it gives the minimum amount of surfactant 
needed to neutralize the DNA charge for the transfection across the cell membrane^^. Since 
the cationic surfactants and lipids are quite toxic, in vivo applications require that the 
amount of amphiphile used in a transfection be as low as possible^. 

The CAC depends on the amount of salt inside the solution and on the specifics of 
the polyelectrolyte and the surfactant to be used. To calculate it, we first obtain the full 
adsorption isotherm, from which the surfoplex formation and the CAC are identified as the 
point of the cooperative binding transition characterized by a sharp rise in the adsorbed 
fraction of the amphiphile. We note that in general this can be either a first or a second 
order phase transition, or simply a sharp crossover. 
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FIG. 2: Density of free surfactant at the transition as a function of salt concentration, log^^Q p'^ (M) x 
Cs (M) for Za = 12, 13, 14, 15. The fihed symbols are the experimental results^.^ and the empty 
symbols are the present theory. The dashed lines are only a guide to the eyes. 

The theory was compared with the experiments of Malovikova et al.^^. who measured 
the adsorption isotherms for sodium dextran sulfate for various concentrations of sodium 
chloride and four different kinds of cationic surfactant: undecyl-, dodecyl-, tridecyl-, and 
tetradeculpyridinium. The four surfactants studied differ by the length of their hydrocarbon 
tail (-2a — 1 = 11, 12, 13 and 14) and, consequently, have different values of the hydrophobicity 
parameter x- Theoretically—, we expect x to be a linear function of Za- The experiments 
were performed at 30°C. The diameter of sodium dextran sulfate is = 7 A and the 
distance between the consecutive charged monomers is 6 = 2.5 A ^ 2.8). The binding 
isotherms obtained experimentally show the formation of surfoplexes when the surfactant 
concentration reaches the CAC, C^{za,Cs,Z,pp), which is strongly dependent on the con- 
centration of all the species and on the type of surfactant used. The CAC was found to 
increase with the salt concentration. The experimental data also showed that the CAC de- 
creased rapidly with the length of the surfactant hydrophobic tail, see Figure 2. To compare 
with the experiment, the theoretical binding isotherms were computed by the procedure de- 
scribed in the previous section. The number density of polyelectrolyte monomers was taken 
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FIG. 3: The hydrophobicity parameter x versus Za obtained by fitting the experimental data of 
Malovikova et al.^. The dashed line is the linear fit. 

to be the same as in the reference^ . The diameter of the polyion counterions, of salt ions, 
and of the surfactant monomers were all taken to be = 9 A. This was done in order to 
simplify the calculations and to minimize the number of adjustable parameters. The final 
results were only weakly dependent on the precise value of Oc. Each polyion was (arbitrarily) 
taken to have Z = 3 000 charged monomers (experiments do not provide the length of the 
polyions, but only the number density of the polyion monomers Zpp). The hydrophobicity 
parameter x "was obtained by performing a fit of the theoretically calculated CACs to the 
experimental data. For each amphiphile type, a unique value of x was used to obtain the 
dependence of the CAC on the concentration of salt. Figure 2 shows that the theory agrees 
quite well with the experimental data of Malovikova et al. for all 4 surfactants. Table 1 
provides the best fit value of the hydrophobicity parameter x, and Figure 3 shows that x is, 
indeed, a linear function of the surfactant length, as was expected theoretically^. The best 
fit of X yields 

X = -1.25^, + 9.5, (11) 

in units of ksT. 

The dependence of the CAC on Cg for various polyion lengths Z at fixed monomer 




FIG. 4: Total density of surfactant at the transition (CAC) as function of salt concentration, 
(M) X Cs (M) foi pm = Zpp = 5x 10-4 M kept fixed and Z = 1 000 (solid line), Z = 2000 
(dashed line), Z = 3000 (dotted line). 

concentration pm — Zpp is illustrated in Figure 4. We see that that, as expected, the CAC 
does not depend individually on Z and pp, but only on the product of the two — the total 
number density of the polyion monomers Pm- Furthermore, from Figures 4 and 5 we see 
that the CAC becomes a linear function of Cs for sufficiently large salt concentrations and 

sufficiently low monomer density, 

CI = 6o {Za. Zpp) + h {Za, Zpp) Cs . (12) 

The amount of salt needed to reach the linear regime depends on the type of surfactant 
used — more salt is needed for surfactants with longer hydrocarbon tails. We now explore 
the dependence of the CAC on the concentration of polyion monomers p^ inside the linear 
regime. Figure 5 show that while the absolute values of the CAC in the linear regime de- 
pend on p„i, the slope is insensitive to the precise concentration of the polyion monomers. 
This implies that bi {za, Zpp) is independent of Zpp so that for a given polyelectrolyte solu- 
tion bi {za,Zpp) — ai{za)- Furthermore, we note that the CAC lines for uniformly spaced 
monomer concentrations are equidistant, see Figure 5. This implies that the coefficient 
bo{za,Zpp) must be a linear function of Zpp, so that bo{za,Zpp) — c{za) + ao{za)Zpp. The 
coefficients 60 ^-nd bi obtained using the list square fits of the linear regime shown in the Ta- 
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FIG. 5: The total density of surfactant at the transition as a function of the salt concentration, 
(M) X Cs (M) for Z = 1 000 kept fixed and = 1 x lO""^ M, 2 x 10"^ M, 3 x lO""^ M, 4 x 10"^ M, 
5 X 10^^ M from bottom to top. The dashed lines are the linear adjust, eqs. ()13p and (jl4p . 

ble 2, clearly support these observations. We thus conclude that for a given polyelectrolyte 
solution, at sufficiently high salt concentrations and sufficiently low monomer densities, the 
CAC is a bilinear function of Z pp and Cs, 

CI = c{za) + ao{za)Zpp + ai(z„) Cs , (13) 

with the coefficients dependent only on the length of surfactant used. This is the fundamental 
result of the present paper. Equation (IT^ should be valid for all rigid polyelectrolyte surfac- 
tant systems at sufficiently high salt concentrations and sufficiently low monomer densities. 
For the specific case of the association of undecyl-, dodecyl, tridecyl-, and tetradeculpyri- 
dinium surfactants with anionic polyelectrolyte sodium dextran sulfate, the coefficients c{za), 
ao{za), and ai{za) can be obtained by interpolating the data presented in the Tables 2. 
We find 

ai (za) = 46 012.89 e-^-^^^" , 
ao{za) = 0.25 + 0.028 
c (za) = -7.47 X 10"^ + 6.57 x 10"^ Za . (14) 
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TABLE II: Values of 6o, h in = 6o + h Cs with Z = 1 000. 

The results for the CAC obtained with the eqs. fll3p and (fT^ are shown as the dashed 
hnes in the Figure 5. 

IV. THE CONCLUSIONS 

We have presented a theory which quantitatively accounts for the cooperative adsorption 
of charged surfactants to rigid polyelectrolytes. The theory is based on the fundamental 
ideas of Debye, Hiickel and Bjerrunti. The only free parameter is the magnitude of the 
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hydrophobic interaction between the condensed surfactants x- agreement with the theo- 
retical expectations, this parameter is found to be a hnear function of the surfactant length 
Za- For high salt concentrations and sufficiently low monomer densities, we find that the 
CAC is a bilinear function of the monomer Zpp and salt Cg concentrations, with the co- 
efficients dependent only on the length of surfactants used. For the specific case of the 
association of the undecyl-, dodecyl, tridecyl-, and tetradeculpyridinium surfactants with 
the anionic polyelectrolyte sodium dextran sulfate, the coefficients of the bilinear form are 
calculated explicitly and presented in eqs. (fT3|l and (fT4|) . 

We note that the present theory does not take into account the non-electrostatic specific 
interactions between the polyions and surfactants. In view of a good agreement between 
the theory and the experiment, however, we conclude that such interactions are not very 
significant, at least for the present polyelectrolyte surfactant system. It would be of great 
interest to test the the predictions of the theory, in particular the bilinear dependence of 
the CAC on the salt and the monomer concentrations, on other polyelectrolyte-surfactant 
systems. 
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